Controlled-NOT logic gate for phase qubits based on conditional spectroscopy 
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A controUed-NOT logic gate based on conditional spectroscopy has been demonstrated recently 
for a pair of superconducting flux qubits [Plantenberg et al, Nature 447, 836 (2007)]. Here we 
study the fidelity of this type of gate applied to a phase qubit coupled to a resonator (or a pair of 
capacitively coupled phase qubits). Our results show that an intrinsic fidelity of more than 99% is 
achievable in 45ns. 

PACS numbers: 03.67.Lx, 85.25.-j 
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I. INTRODUCTION 

Approaches to the development of a large-scale quan- 
tum computer face numerous practical challenges 
0|. One such challenge, the construction of a robust 
controlled-NOT (CNOT) logic gate, has been explored 
from several perspectives [4-10]. In this work we analyze 
a CNOT gate based on conditional spectroscopy for a 
superconducting phase qubit coupled to a resonator. A 
related construction has already been demonstrated for a 
pair of superconducting flux qubits and a similar con- 
cept has been explored in the context of NMR quantum 
computing Although this approach can be applied 
to a variety of physical systems, the fidelity observed in 
the experiment of Ref. |llj] is not sufficient for practical 
use. Here we propose a method to improve the intrin- 
sic fidelity and we then calculate the optimal fidelity as 
a function of total gate time. The phase qubit coupled 
to resonator system we consider is relevant to the UCSB 
Rezqu architecture, [isj 

The idea behind a spectroscopic CNOT gate is simple 
and has a wide range of applicability: A tt pulse is ap- 
plied to the target qubit with a carefully selected carrier 
frequency Wrf . The carrier frequency is close to the qubit 
transition frequency given that the attached control res- 
onator is in the "on" or |1) state, which in the \qr) basis 
is 
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The Rabi frequency has to be smaller than the detuning 
to the "off transition at 
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A direct ^ coupling between the devices would 
of course generate a difference in Won and Woff, but in 
the phase qubit plus resonator system — which has no di- 
rect (7^ (8) cr^ coupling — such an interaction is generated 
by level repulsion from the noncomputational |2) states. 



The difference |a;on — Woff | characterizes the sensitivity of 
the conditioning effect and determines the speed of the 
resulting gate. 

When the qubit and resonator are detuned by an 
amount larger than the coupling between them, they be- 
come weakly coupled. In this limit, the sensitivity for 
current devices is limited to a few MHz, which is not suf- 
ficient for practical application. Therefore to amplify the 
sensitivity we adiabatically bring the target qubit to a 
suitable point near resonance with the control resonator, 
and drive the qubit while it is strongly coupled with the 
resonator. After performing a tt pulse the qubit is adia- 
batically detuned from the resonator. 

Two main sources of intrinsic errors exist in this ap- 
proach: Although we set the carrier frequency Wrf to a 
value such that we have a tt pulse in the qubit when the 
resonator is in there is a small probability for the 
qubit to get rotated even if the resonator is in |0). The 
second error comes from the fact that, since we are driv- 
ing the qubit, it is possible to have leakage to the qubit 
1 2) state. The fidelity will reach its maximum value when 
both of these errors are minimized simultaneously. We 
use the DRAG method [l4] to suppress the error due to 
leakage and adjust all other parameters by optimization. 



II. CNOT DESIGN 



Hamiltonian 



In the basis of uncoupled qubits, the Hamiltonian of a 
qubit capacitively coupled to a resonator (assuming har- 
monic eigenfunction of 2-states) is given by (suppressing 

h), 
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H — Hq + iJrf + Hii 



(3) 
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where, 
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where, e, w and A are qubit frequency, resonator fre- 
quency, and anharmonicity of the qubit, respectively. 
Q., Wif and 4> are Rabi frequency, carrier frequency, 
and phase of the microwave pulse, g is the (time- 
independent) interaction strength between qubit and res- 
onator. 
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FIG. 1: (Color online) A plot of first six energy levels of 
a phase qubit coupled to a resonator a.i g/h — lib MHz. 

The first six energy levels of the system (obtained 
numerically) are shown in Fig.(IT]). We use coupling 
f;/ft,= 115 MHz, resonator frequency uj/h — 6.5 GHz and 
anharmonicity of the qubit A/h = 200 MHz. The simula- 
tions discussed below are carried out in a frame rotating 
with the instantaneous frequency of the qubit. 



B. CNOT Protocol 

Eigenstates of the full Hamiltonian reduce to the eigen- 
states of the uncoupled Hamiltonian far away from the 
resonance. We denote the first six eigenstates of the full 
Hamiltonian by Eqq, Eqi, Eiq, Eq2, Eh, £'20 and define a 



conditional control sensitivity Sc and leakage sensitivity 
Si as 



Sc — \{Eii — Eqi) — {Eio — Eoo)\, 
Si = \{E2i-En)-{En-Eoi)\. 



(5) 



The conditional control sensitivity is the (magnitude of 
the) difference between Won and lOqS ■ Leakage sensitivity 
is the anharmonicity of the target qubit when control 
resonator is on. In order to achieve a high fidelity both 
of these quantities need to be maximized (by varying the 
detuning e — cu). 
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FIG. 2: (Color online) A plot of control and leakage 
sensitivity vs. qubit frequency for g/h = 115 MHz. 



A plot of these sensitivities is shown in Fig. [5] Peaks 
in the control sensitivity (resulting from expected anti- 
crossings) at detuning equal to zero and A conspire to 
give the maximum at 100 MHz detuning shown in the 
blue curve of Fig. [51 Operating near 100 MHz detun- 
ing, however, leads to poor performance because of the 
large leakage error there. Better operation points exist 
near -100 and 200 MHz detuning; we shall make use of 
the latter. Fig. ([3]) shows the behavior of sensitivities vs. 
coupling at {e — ijj)/h — 215 MHz. 

To implement a CNOT gate we begin with a strongly 
detuned qubit-resonator system. Then the qubit is adia- 
batically tuned into resonance with the resonator, driven 
with a TT pulse, and finally detuned. In the \qr) basis this 
protocol ideally produces 
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= SWAP X CNOT X SWAP, 



(6) 
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FIG. 3: (Color online) A plot of control and leakage 
sensitivity vs. coupling. 



where 



SWAP = 



10 
10 
10 
1 



(7) 



is the swap gate. To obtain this target we also perform z 
rotations on the qubit before and after the sequence de- 
scribed above, with angles determined by optimizations. 



C. Local Rotations with DRAG 

As far as leakage outside the computational subspace 
is concerned, we can consider Eoi,Eii and E21 to be a 
single 3-level quantum system where we are interested in 
local rotations between the first two levels and therefore, 
in order to suppress the errors due to leakage to the third 
level, local rotations are performed with a DRAG pulse 
[14], up to 5"^ order. In order to do a local rotation of 
angle 9 about the x axis, we set the Rabi pulse of the 
Hamiltonian to be 



r!(i)cos(0(<)) = /e 
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where, A = for our case and the first two equations 
give the amplitudes of cos(a;rft) and sin(a;rft) quadra- 
tures. Here uJc is found via optimization around I-E01--E11 1 
and feit) is chosen to be a Gaussian function (vertically 



shifted) such that fg{Q) = fe{tg) = and j fe{t)dt = 6*, 



tg being the time required to perform the local rotation. 

III. GATE OPTIMIZATION 

In this section we describe the sources of error and our 
methodology to compute the intrinsic fidelity curve. As 
mentioned above, there are two major intrinsic sources 
of error. In our simulation we use the full nine-level 
Hamiltonian to compute the fidelity curve and use DRAG 
to suppress the leakage error while the control error is 
treated with an adjustment of other controllable param- 
eters via optimization. 

The state-dependent process fidelity between two 
quantum gates U and f/targot is given by. 



F(C/targct,C/,|V)) 
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The fidelity measure we use in this work is the process 
fidelity averaged over the four dimensional Hilbert space, 
that can be shown [l^ [l^ to be exactly equal to a trace 
formula given by. 



i^(f/target,f/) 



Tr((7;7t)+ Tr (C/tc/) 
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(10) 



where [/target is always four dimensional for a two-qubit 
operation and U is the time evolution operator projected 
into the four dimensional computational subspace of the 
entire Hilbert space. This definition assumes that C/targct 
is unitary, but U does not have to be (because the fi- 
nal evolution operator after projection into the four- 
dimensional computational subspace is not necessarily 
unitary). 

In order to compute the fidelity curve, we perform a 
8 dimensional optimization over the control parameters 
for a given total gate time and obtain the best fidelity. 
The 8 control parameters optimized here are coupling (g), 
qubit frequency in the intermediate state (see FigH]), car- 
rier frequency of the Rabi pulse (wc), angle of rotation by 
the Rabi pulse, duration of each ramp pulse (iramp) and 
three other parameters related to the areas of the gaus- 
sian envelope of DRAG pulse. The result obtained from 
such an optimization corresponds to a single point in the 
plot of Fig.([3]). In order to help to control the phases 
of the matrix elements developed in the time-evolution 
operator, we also attach a pre and a post z-rotation of 
the qubit. We observe from Fig. ([2]) and Fig. ([3]) that a 
good range for driving point (e — cj) should be ~ 200-250 
MHz and a good range for coupling strength should be 
~ 100-125 MHz. These numbers are used as guess values 
of our multidimensional optimization search. 

Fig.([3]) shows the fidelity curve obtained from op- 
timization and Table jll shows corresponding coupling 
strengths, duration of each ramp and driving points. Our 
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TABLE I: Optimum values of parameters for CNOT 
gate. 



^gatc C^^O 


tram (uS ) 


q (MHz) 


e-uj (MHz) 


Fidelity [%] 


33 


7.1 


118.5966 


235.1804 


92.4021 


35 


7.3 


107.5055 


218.8692 


97.2223 


37 


7.0 


123.8367 


237.4236 


96.4515 


39 


7.3 


113.2576 


217.7214 


98.7965 


41 


7.4 


107.7414 


218.2388 


98.5218 


43 


7.2 


117.6687 


218.6669 


98.5886 


45 


7.5 


107.3625 


199.2662 


99.3914 


47 


7.5 


104.0402 


204.6788 


99.2326 


49 


7.3 


113.8836 


208.3964 


99.1864 


51 


7.5 


107.6466 


202.2407 


99.3836 



result shows that the intrinsic fidelity can be pushed to 
99% within 45 ns gate time. The remaining error comes 
from the adiabaticity of the ramp pulses and leakage 
outside the computational basis states, for example be- 
tween 1 11) and 1 20) (in |qr) basis) at driving point where 
e « w + A. 
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FIG. 4: (Color online) A plot of fidelity vs. total gate 
time in nanoseconds. Red points are obtained via 
optimization and the blue line is an interpolation. 



As an example, we show the change of qubit frequency 
(in GHz) over time (in nanoseconds) in Fig. ([5]) for the 
CNOT having total gate time — 45 ns. while the res- 
onator frequency is always fixed at 6.5 GHz and pre and 
post z-rotation angles are found to be impost = —1.0915 
and i?prc — 0.5442 radian for this case. We use lin- 
ear pulse for ramps and a gaussian envelope is used for 
DRAG. 
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FIG. 5: (Color online) A plot of qubit frequency vs. 
time for CNOT at total gate time=:45 ns. 



IV. CONCLUSIONS 

We have computed intrinsic fidelity of a CNOT gate 
based on conditional spectroscopy approach and have 
shown that it is possible to achieve greater than 99% fi- 
delity within an experimentally practical time scale of 
45ns. However, this design req uires a large coupling 
strength, so tunable coupling |17l| would probably be re- 
quired in a multi-qubit system. Although our analysis 
assumed a phase qubit and resonator, the design also ap- 
ples to capacitively coupled phase qubits, where a slightly 
higher fidelity would be expected because of the addi- 
tional anharmonicity. 
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